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(54) Method of estimating time offsets in a position determining system 



(57) The invention relates to a position determining 
system, for receiving digital telephone signals transmit- 
ted by a number of transmission sources BTS. In par- 
ticular there is described a method of estimating the off- 
set in time of receipt of a broadcast signal received at 
two locations (CBU,CRU), where the signal received at 
the one location (CRU) may have been subject to cor- 
ruption by multipath effects. To overcome this problem, 
the signals received at the one location are auto-corre- 
lated and the signals received at the other location are 



auto-correlated. The signals received at the two loca- 
tions are cross-correlated and a template comprising 
that portion of the auto-correlation of the signals re- 
ceived at the other location corresponding to the nega- 
tive time axis, and that portion of the auto-correlation of 
the signals received at the one location corresponding 
to the positive time axis is constructed. The offset at 
which the template best fits the measured cross-corre- 
lation of the signals received at the two locations is 
measured as an estimate of the time offset between the 
signals received at the two locations. 
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Description 

[0001] This invention relates to a method of estimating time offsets in a position determining system using radio or 
other broadcast transmissions. More particularly, it relates to a system which uses the Global System for Mobile Com- 

5 munications (GSM) or other digital transmission systems. 

[0002] In our European Patent EP-B-0 303 371 we describe a radio navigation and tracking system, now known as 
'CURSOR', which uses the spatial coherence of the signals from several radio transmitters to determine the position 
of a roving receiver. Its principles are explained in the above patent specification where it is shown how the signals 
received directly by the roving receiver (rover) are compared with those received by a fixed base station whose location 

10 is known (base) to determine their phase difference, and hence the difference in range of the base and rover from each 
transmitter. Three such measurements made on independent transmitters are needed for navigation and tracking in 
two dimensions to fix the position of the rover relative to the base station and network of transmitters. The unknown 
quantities calculated for each new position are the spatial x and y co-ordinates of the rover together with the phase 
offset between the local oscillators in the equipment of the two receivers. A further patent specification, WO 94/28432, 

*5 shows how the same principles may be applied in tunnels and other shielded spaces such as underground car parks. 
Duffett-Smith and Woan (Journal of Navigation, 45, 157, 1992) describe a particular implementation of the phase- 
measuring system in which the signals from three or more medium-wave AM public broadcast stations are used to 
track the position of a vehicle in and around Cambridge, UK, at speeds up to 110 km h" 1 , with a precision of about 5 
m. One of the advantages of such a system lies in the fact that no additional costly infrastructure of coherent radio 

20 transmitters needs to be set up for CURSOR operation. On the contrary, CURSOR is able to use the signals from any 
independent radio transmitters set up for any purpose. 

[0003] EP-B-0 303 371 also explains how the wider bandwidth signals of modulated transmissions may be used to 
measure the time difference between the signals received from each transmitter at the base station and the rover. In 
this case, the position of the peak in the cross correlation can be used as an estimator of the time difference between 
25 the two received signals, and hence the difference in distance from the transmitter of the base and rover. As with the 
phase-measuring system, three such measurements made on three widely-spaced transmitters suffice to calculate 
the spatial x and y co-ordinates of the rover together with the time offset between the oscillators in the two receiving 
stations. 

[0004] Mobile telephones are increasingly using GSM and other digital techniques, and it would be an advantage to 
30 add time-measuring CURSOR positioning technology to provide additional services for users. However, signals radi- 
ated by digital telephone transmitters are complex and therefore there are problems to be overcome in attempting to 
combine the technologies. 

[0005] According to the parent application (EP-A-088071 2) of this application there is provided a position determining 
system, for receiving broadband signals transmitted by a number of transmission sources (BTS) equal at least to the 
35 number of dimensions in which the movement of a roving object (R) is to be monitored, the system comprising a pair 
of receiving stations, in use the first of the receiving stations (CBU) being at a known position and the second (CRU) 
being located on the roving object; 
characterised in that 

each of the receiving stations (CBU, CRU) is arranged to receive the signals from the respective transmission 
40 sources (BTS) substantially simultaneously, and record a portion of each of the received signals, and each receiving 
station including means for transmitting a link signal (L v Lg) containing information about the recorded portions of the 
signals received from the transmission sources; 

and by 

a position determining processor (CPP) arranged to receive the link signals (L 1( L 2 ) and compare the information 
45 received from the one receiving station about the recorded portions of the signals received at the one receiving station 
from the transmission sources with the information received at the other receiving station about the recorded portions 
of the signals received at the other receiving station from the transmission sources, and to determine the time delay 
between the respective signal portions received at both receiving stations in order to determine the position of the 
roving object (R). 

so [0006] The receiving stations may receive the signals from the respective transmission sources sequentially and in 
the same sequence as each other. 

[0007] The parent application also includes a method of determining the position of a roving object (R), the method 
comprising transmitting broadband signals from a number of transmission sources (BTS) equal at least to the number 
of dimensions in which the movement of the roving object is to be monitored; receiving the signals at a pair of receiving 
55 stations (CBU.CRU), in use the first (CBU) of the receiving stations being at a known position and the second (CRU) 
being located on the roving object (R), 
characterised in that 

each of the receiving stations (CBU, CRU) receives the signals from the respective transmission sources (BTS) 
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substantially simultaneously, records a portion of each of the received signals, and passes a link signal (L 1t L 2 ) from 
to a position determining processor (CPP), the link signals containing information about the recorded portions of the 
signals received at the respective receiving station from the transmission sources; and 

the position determining processor (CPP) compares the information received from the one receiving station (CRU) 

5 about the recorded portions of the signals received at the one receiving station from the transmission sources (BTS) 
with the information received from the other receiving station (CBU) about the recorded portions of the signals received 
at the other receiving station from the transmission sources, and determines the time delay between the respective 
signals received at the receiving stations in order to determine the position of the roving object. 
[0008] According to the present invention there is provided a method of estimating the offset in time of receipt of a 

10 broadcast signal received at two locations, where the signal received at the one location may have been subject to 
corruption by multipath effects, the method comprising auto-correlating the signals received at the one location; auto- 
correlating the signals received at the other location; cross-correlating the signals received at the one and other loca- 
tions; constructing a template comprising that portion of the auto-correlation of the signals received at the other location 
corresponding to the negative time axis, and that portion of the auto-correlation of the signals received at the one 

15 location corresponding to the positive time axis; and measuring the offset at which the template best fits the measured 
cross-correlation of the signals received at the one and other locations as an estimate of the time offset between the 
signals received at the two locations. 

[0009] In some systems, each of the receiving stations is also arranged to receive a second signal from one or more 
of the transmitters, said second signal being used to enable compensation for changes in equipment offsets during 
20 receipt of the sequentially received signals. 

[0010] The position determining processor may be co-located with one of the receiving stations or may be located 
remotely. 

[0011] In one embodiment, the link signal from one of the receiving stations is passed to the other receiving station 
and from the other receiving station to the position determining processor. 
25 [0012] Preferably, for locating purposes, a signal providing information about the position of the roving object is 
passed from the position determining processor to at least one of the receiving stations. 

[0013] The system may further comprise one or more monitoring stations and a signal providing information about 
the position of the roving object may then be passed from the position determining processor to one or more of the 
monitoring stations. 

30 [0014] To provide enhanced functionality, a database server may be connected to the position determining processor, 
the database server containing data elements relating to a plurality of known positions; and the system further com- 
prises means for passing information about the position of the roving object determined by the position determining 
processor to the database server; means for retrieving data elements related to the position determined by the position 
determining processor; and means for passing the data elements to one of the receiving stations or to one or more of 

35 the monitoring stations. 

[0015] The or each receiving station or the or each monitoring station preferably includes a display and the position 
of the roving object is displayed on the display which may be a dot matrix display. 

[0016] The database server may contain graphical information and the graphical information is passed to a receiving 
station or monitoring station and displayed on the display to indicate the position of the roving object. 

40 [0017] The transmission sources and receiving stations preferably comprise components of a digital cellular tele* 
phone network such as a GSM network. It is advantageous if the receiving stations monitor the strength of the signals 
of plural transmission sources and select a plurality of sufficiently strong signals for reception. 
[0018] The quasi-synchronisation between the signals received from the respective transmission sources is prefer- 
ably achieved by monitoring receipt of a specific portion of the transmitted signals, but may be achieved, alternatively, 

45 by means independent of the transmitted signals such as a local clock signal. 

[0019] The system preferably includes a series of regionally, nationally, or even internationally networked fixed re- 
ceiving stations to provide wide coverage and usage of the positioning system. 

[0020] The receiving stations may be capable of receiving two or more channels simultaneously and it may be ad- 
vantageous to repeat the reception of signals from plural transmission sources. 

so [0021] The description below sets out how the principles of the time-measuring CURSOR system may be applied 
to a digital radio network, such as the GSM telephone system, to enable the position of a receiver, such as a portable 
telephone handset, to be measured relative to the transmitter network. The equipment used in such a network already 
incorporates most of that required for position fixing operation; so that its implementation may be achieved with little 
extra cost. The precision of each position measurement is proportional (amongst other factors) to the inverse of the 

55 bandwidth of the signals; for a GSM single channel of 200 KHz bandwidth, the precision is about 50 m. Some improve- 
ment on this figure may be obtained in cases where more than three GSM transmitters can be received at both base 
and rover; they may all then be used in the position determining process. 

[0022] The principles of operation of a position determining system according to the present invention, and a specific 
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application to GSM technology, will now be described with reference to the accompanying drawings and included 
tables, in which:- 

Figure 1 is a diagram of the basic system elements; 
s Figure 2 is a diagram used to define the system co-ordinates; 

Figure 3 illustrates calculation of cross correlation functions; 

Figure 4 illustrates a GSM network incorporating the time-measuring technology of the invention, and identifies 
the logical elements thereof; 

Figure 5 illustrates the basic circuit components of a GSM handset; 
10 Figure 6 illustrates a particular method of overcoming the effect of multipath propagation on measuring relative 

time delays; 

Figure 7 (a, b, c) shows various different configurations of the elements in a GSM system incorporating the inven- 
tion; 

Figure 8 illustrates how position information may be displayed graphically on a small liquid crystal display; 
15 Figure 9 illustrates the arrangement of elements and test sites in a trial carried out in Cambridge, UK; 

Table 1 shows how the intrinsic measurement error varies with the signal to noise power ratio and the number of 
bits recorded; and 

Table 2 gives some results from tests of a position determining system carried out in Cambridge, UK. 



20 Principles of operation 



[0023] Let us suppose that a transmitter, A, is broadcasting modulated radio signals within a band of frequencies of 
width Av centred on frequency v 0 . The signals are received by a fixed base station, O, and a roving receiver, R, as 
shown in Figure 1 . We define Cartesian co-ordinates (x, y) with respect to axes centred on the base station, as shown 
25 in Figure 2. The axes can have any orientation, but might conveniently be set so that the y-axis lies along the local 
mapping grid north-south direction, a - {a„ a y ) is a vector which defines the position of the transmitter, and r = (x, y) 
defines the position of the rover. The vector r - a = (x-a x , y-ay) completes the triangle OAR. 

[0024] Transmitter A transmits signals continuously. At a prearranged moment the radio receivers in O and R begin 
receiving and recording a short length of the transmissions. This moment may for example be signalled by the arrival 

30 of a trigger transmitted from A, or a trigger derived from the normal signal traffic radiated by A, or a trigger generated 
locally. The last of these may be achieved, for example, by using a pulse of short duration derived from high-precision 
clocks at R and O which have been synchronised earlier. In a digital version of this application, the signals might first 
be converted from V Q to baseband, digitised, and then recorded in dynamic memory. Although O and R record as near 
simultaneously as possible, there will nevertheless be a time offset between the two records which is caused both by 

35 the difference in the distances AO and AR, and by the synchronisation error between the arrival of the trigger at O and 
R. Let At be this time offset; it is given by 



40 



At = 



[(* 



a/ + y- 



2.% 



+ £ 



where e is the synchronisation error and t) is the speed of the radio waves. 

[0025] An estimate of At may be obtained from the cross correlation of the recordings made at O and R. In Figure 1 
there are shown links, L1 and L2, which represent means of conveying representations of the recordings made at each 
receiving station to a position determining processor CPP. The nature of the links, whether they convey information in 
near real time or impose delays, is of no consequence in principle. All that matters is that duplicates of the recordings 
made by each receiving station are assembled together at the position determining processor for comparison with 
each other. 

[0026] The cross correlation of the two recordings is now performed by the position determining processor using any 
convenient means. In the case of digital signals, this may best be carried out by a microprocessor which calculates 
estimates of the function at discrete intervals of time offset, x, corresponding to the sampling interval, f s , as shown by 
the dots in Figure 3. It is unlikely that the peak of the function corresponds with a particular sample, so the microproc- 
essor must also fit an interpolation function, shown by the dashed curve, in order to obtain the best estimate of the 
position of the peak. The value of r corresponding to the peak is an estimator of At. 

[0027] This process of quasi-simultaneous recording, conveyance over one or more links, assembly of copies at one 
place, and estimation of the value of At from the cross correlation, must be performed on at least two or three spatially 
distinct transmitters. If there are means at O and R of synchronising sufficiently accurately, then just two measurements 
would generally suffice; otherwise, if £ is unknown as is generally the case, at least three measurements are needed. 
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10 



15 



55 



Let A, B, and C be the three transmitters in the latter case, at vector positions a = (a x , ay), b - (b x , by), and c = (c^, Cy). 
Then the three measurements of Af are given by 

[(x-a/+y-a y f?-(a 2 x + a 2 / 
AL = - — + e 



Af l(x - bf + (y - b y ff - (b\ + b 2 / | ^ 



and 



Atc Jx-c/+(y-c y f?-(4 + c 2 / + e (1) 



[0028] These three non-linear equations can be solved to evaluate x, y, and e, and hence the position of the rover 
is determined. Actually, the solution of these three equations is generally ambiguous with two possible pairs of (x, y). 
20 Often these are sufficiently far apart that it is obvious which is the correct result (for example, only one might lie within 
the triangle formed by the three transmitters), but if the ambiguity must be resolved automatically then an extra trans- 
mitter must be monitored. 

[0029] The accuracy of the method depends on the error in the estimator of Af derived from the cross correlation, 
and there are three main factors. First, the width of the cross correlation affects the accuracy of the result, since a 

25 wider function provides a less sharply-defined peak. The width is inversely proportional to the bandwidth, Av, of the 
transmitted signals. Hence, wider bandwidth transmissions provide for more accurate intrinsic position determination. 
It should be noted here that the "bandwidth" could encompass all the signals radiated from a particular radio mast and 
not just one particular channel. For example, if ten adjacent channels of 200 KHz width each were active within a total 
bandwidth of 2 MHz, it may be possible, using a wide-band receiver or one which can receive simultaneously on more 

30 than one channel, to make use of all ten at once and thereby obtain the accuracy corresponding to a 2 MHz bandwidth. 
If some of the channels were inactive, it might nevertheless still be possible to obtain the full accuracy corresponding 
to the separation of the two channels furthest apart. 

[0030] Secondly, the signal-to-noise power ratio, r, of the signals is important. Let us assume that we are recording 
one-bit samples of the received signals converted first to baseband. Each sample is a 'one' or a 'zero' depending on 
35 whether the signal is greater than or less than zero at the instant of sampling. We also assume that the signals are 
contained wholly within a uniform band of frequencies extending from zero to Av, and that they are sampled at the 
Nyquist rate equal to 2Av. The error in the estimation of the position of the peak, Ax, is related to the average number, 
q(r), of consecutive samples which contains just one error by the approximate expression 

40 



(q(r)-2)AvjN 

where N is the total number of samples used in the cross correlation. The minimum value of q is 2 since even when 
45 there is no signal, the one-bit samples have an equal chance of being right as wrong. Here we have assumed that the 
cross correlation function is triangular in shape rather than the sinc-squared function expected for noise signals and 
square band-pass filters. An error of Ax adds in quadratically as a position error, Ax, where Ax= v Ax. The relationship 
between q(r) and rcan be obtained from the following expressions on the assumption that the signals are noise-like: 

50 i 
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With Av set to 200 KHz, the results shown in Table 1 are obtained. 



Table 1 
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[0031] As noted above, the errors on the measurements from each transmitter,say Ax 1( Ax 2 , Ax 3 , must be added in 
quadrature with each other and with any other errors in the estimation of the position of the rover. It is interesting to 
note that quite good results can be obtained even in conditions of poor reception. 

[0032] Thirdly, and usually the most important in practice, there is the error incurred by not having an accurate knowl- 
edge of the paths by which the signals reach the rover. Multipath propagation spreads out the cross correlation, making 
it harder to estimate the position of the peak. It may also result in a multi-peaked cross correlation with the desired 
peak having a lower amplitude than others. If all the signals arrive by indirect routes, there may be no peak at ail 
corresponding to the line-of-sight propagation path. It should be noted, however, that multipath propagation always 
results in a delay of the signals compared to the direct path. Provided that the base station antenna is in the clear 
above the surrounding clutter so that it receives the most direct signals only, then delayed signals at the rover always 
appear to the later side of the peak of the cross correlation. (In these circumstances it is possible to alleviate the effects 
of the multipath propagation as explained later) Having fitted an appropriate interpolation function to the discrete 
samples, it is important to choose the smallest value of r at which there is a significant amount of signal as that value 
of Af to be used in the calculation of the position, rather than the position of the peak itself. 

Application to GSM 

[0033] The signals radiated by the transmitters in a GSM network are complex. The flexibility and capacity built into 
its design are such as to render it difficult, if not impossible, to predict precisely how the allocated spectrum (near 900 
MHz and 1 800 MHz in Europe) will be used at any moment. The frequency band is split into a number of radio frequency 
(RF) channels of 200 KHz width, each of which carries frequency modulated (FM) signals divided in time into a sequence 
of frames. The basic unit is the so-called time division multiple access (TDM A) frame lasting 4.615 ms and further 
divided into 8 time slots. Each time slot carries 156.25 bits at a rate of 270 Kbits s _1 and may represent a 'normal burst' 
of data and training bits, a 'frequency correction burst* of fixed pattern, a 'synchronisation burst' of data and synchro- 
nisation bits, or an 'access burst' with a synchronisation sequence and data. Each of these bursts also carries header, 
tail, and guard bits. How many of the time slots are being used at any moment in a given frame, and how many of the 
RF carriers are being broadcast from a given transmitter, depends on the way the system has been set up and on the 
amount of traffic at that moment. However, in the minimum case when all is quiet one of the RF carriers will always be 
alive, carrying the so-called broadcast control channel (BCCH; a logical channel), polling for handsets in its cell by 
broadcasting one access burst in every TDMA frame. We can therefore rely on there being at least one FM broadcast 
from each transmitter mast with a bandwidth of about 200 KHz, and this we can use for position location. 
[0034] A regional GSM network incorporating a positioning system according to the present invention might be as 
shown in Figure 4. Here we see base receivers, hereinafter also referred to by the acronym CBU, labelled CBU A , 
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CBU Bt CBU C etc., at every GSM transmitter A, B, C etc., able to receive not only the signals from their own local 
transmitters but also those from at least two other distant transmitters. We also see a number of handsets incorporating 
roving receivers, hereinafter also referred to using the acronym CRU, labelled CRU V CRU 2 etc., active in the region. 
It is important that these handsets can also receive the signals from the same set of distant transmitters as well as 

5 from the local transmitter. As noted above, it is not necessary to achieve a high signal-to-noise ratio for reception of 
the distant signals, but nevertheless the need to receive distant transmitters may limit the usefulness of this system in 
rural areas where cells are very widely separated. Also shown in Figure 4 are the position determining processor (CPP) 
and a service provider database (SPD). This is a device which provides position- related data specific to a request by 
the user of the position determining system. For example, the user may need directions to find the railway station in 

10 an unfamiliar city. The position determining processor would compute the user's position and pass it to the service 
provider database together with the user's request. The database would respond with the required list of instructions. 
The configuration shown here, in which the position determining processor is a stand-alone remote unit, is just one of 
several possible configurations (see below). For example, the CPP and CRU might be combined together inside the 
handset so that position processing is carried out by the handset's own computer. 

is [0035] There needs to be established a trigger code which is broadcast periodically from each transmitter, probably 
within the BCCH logical channel. As noted above, this coutd be a special code set up within the existing GSM framework, 
or it might be a repeating element of the ordinary GSM signals themselves, such as the arrival of a frame number 
ending in three zeroes which occurs every 4 seconds. The arrival of the trigger code causes an active roving receiver 
(CRU) to begin the process of recording the signals from at least three transmitters. There may sometimes be a com- 

20 mercial advantage in making the trigger code unique to a particular operator, or even to a particular handset, thus 
providing one means of charging users for the position determining service. (There are of course many other possible 
ways of charging.) The trigger code also activates the recording mechanism in the base unit (CBU) attached to the 
GSM transmitter. Since most GSM handsets can receive only one RF channel at once, the order in which to record 
the distant signals must first have been established, for example by using cell broadcast or the short message service. 

25 We see below that it might also be necessary to record the local transmitter signals a second time in order to take 
account of drifts over the recording period. If 2048 samples are recorded for each of the two distant transmitters (1 024 
each of I and Q quadrature-phase signals), and 4096 samples for the local transmitter (two sets of 2048), the whole 
process can be completed within a few hundred ms including the time taken to switch between the channels and settle 
on each new frequency, The recordings made by the roving receiver(CRU) and base receiver (CBU) are then sent by 

30 any convenient means to the position determining processor (CPP). The roving receiver could use, for example, slots 
in the GSM signals, or it could set up a data transfer call. The base receiver would probably send its recordings by 
landiineto a remote position determining processor. The processor receives and assembles in memory its own copies 
of the signals received by the base and roving receivers for correlation and position processing. It might also use other 
parameters of the signals recorded by the two receiving stations such as signal strengths. 

35 [0036] Once the data have been assembled, the processor then performs the cross correlation analysis. Standard 
methods may be used, but a preferred method is described later which reduces the effects of multipath propagation. 
Having made estimates of Af a , A^ and A^., the processor solves the equations to evaluate x and y. However, the need 
to receive the signals from the three transmitters sequentially introduces a complication in that the synchronisation 
error, e, is changing all the time and is unlikely to be constant over the period of recording. We can model this error 

40 over short periods as an offset plus a linear slope, i.e.: 

e = e 0 + 

45 where £q and e t are constants, and f is the time. We have a further unknown quantity to evaluate, e 1t and we can do 
this by recording the signals from the local transmitter again for a second time (at both base and rover). We now obtain 
four values of At made at sequential times f 1( t 2 , t z , and f 4 , and equations (1) become: 

so A#/# . [(x~af+(y-a/F-(4+#* 

A ,„, [(X't>/+(y-b y ) 2 F-(t> 2 x + b 2 / 
55 A WW = ~ — +e o + e 1 f 2 
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Af c (f 3 ) = 



l(x-c x ) +(y-c. 



ft-. 



cUc 2 f 



+ 60 + 8^3 



5 and 



10 



The first and last of these equations can be subtracted to give 



15 

from which e-j, can be found. Hence with four estimates of Af equations (2) can be evaluated to recover x and y. 
[0037] Handsets containing multi-channel receivers able to receive all three channels simultaneously do not need 
to repeat the recording of the signals from transmitter A since the rate of drift between the clocks (measured by is 
then of no consequence and does not need to be determined. High-quality single-channel handsets may also be able 

20 to dispense with the repeated recording if the drift in the clock error is small enough to ignore. As noted above, almost 
all the hardware elements needed inside a GSM handset for the purposes of the present invention exist already for 
normal GSM operation. Figure 5 illustrates a typical implementation. The RF signals in the region of 900 MHz or 1 800 
MHz are received by antenna 1 , amplified in the RF amplifier 2, and filtered in the RF bandpass filter 3 ( before being 
mixed to the intermediate frequency, IF, in mixer 4 with a local oscillator signal L0 1( generated by oscillator 5. L0 1 

25 already incorporates the necessary channel switching ability. It is important for the present invention that its frequency 
is phase-locked to that of a master oscillator, MO, usually a crystal oscillator 6. This is also a requirement for normal 
GSM operation so it imposes no additional costs. The IF signals are amplified in amplifier 7, filtered by the bandpass 
filter 8, and converted to baseband in quadrature mixers 9 and 1 0 using the signal, L0 2 , from a second local oscillator 
11 . As with L0 1( this second oscillator also needs to be phase-locked to the master oscillator. However, fine tuning 

30 (automatic frequency control, AFC) is often also incorporated into the loop which must be held constant during the 
recording process. Quadrature outputs are filtered by the lowpass filters 12 and 13, digitised in analogue-to-digital 
converters 14 and 15, and the bit streams labelled I and Q are produced. Position determination then just requires the 
microprocessor 1 6 to make recordings of the bit streams in dynamic ram 1 7 u nder the control of the program 1 8. There 
might be a need for some additional memory. The main modification to enable position determination is in the firmware 

35 program running on the microprocessor, and it might be necessary to upgrade the processor itself to a faster model. 
Similar equipment is needed in the base receiving station. 

[0038] A commercially limiting element in a system according to the present invention is likely to be the means by 
which the signals recorded at the base and rover are transferred to the position determining processor for cross cor- 
relation. It is therefore important to reduce the number of bits transferred as much as possible. As noted above, one- 

40 bit sampling produces adequate results even under conditions of poor signal reception, and although slightly better 
accuracy can be obtained using, say, two-bit sampling, the increase in the link load is likely to be unattractive. Indeed, 
it is probably better to double the number of one-bit samples, and hence double the length of the sampling interval, 
than to use two-bit samples. It may also be necessary to measure both I and Q channels and to perform a complex 
cross correlation, since the base and rover signals are unlikely to be in phase with each other, and the crystal oscillators 

45 are not locked together. Let the quadrature bit streams at the rover (CRU) be represented by / and Q (see Figure 5), 
where 



/= vfjcosfl,; 

50 

and 



V 2 (t)sm(B 1 ) 

55 

where 9 1 is the phase of the signals, and l/ f (f), V 2 {t) are their amplitudes. The corresponding quantities recorded at 
the base (CBU) are 
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and 

5 

o'= i4(f)sinre 2 ; 

The microprocessor must then compute four cross correlations corresponding to //', QQ, IQ, and I'Q, and then seek 
10 the peak in the quantity p, where 



p^Jfll'+QQ'f+OQ-l'Qf 
is The phase difference, 0= QyB 2 , can also be obtained from 

0=tan { TFTQQ i 

20 and this can be used to take out the small frequency offset between the local oscillators in the base and rover so that 
longer integration is possible. 

[0039] Equations (1) and (2) above show how estimates of the time delays, Af, may be used to deduce the position 
(x, y) of the rover and the synchronisation error, e, between the base and rover. Central to this process is the accuracy 
with which estimates of At can be obtained from the cross correlations. As noted above, the major limitation to the 

25 accuracy is likely to be caused by multipath propagation of the signals from the GSM transmitter to the rover. It is 
assumed that the antenna at the base has been installed high enough above surrounding clutterto ensure that multipath 
propagation is small enough to be ignored. In these circumstances, multipath propagation to the rover causes the cross 
correlation profile to be spread out towards longer delay times, since the most-direct signals arrive first. The preferred 
method for obtaining an estimate of At is based on using auto-correlations (ACF) of the base and rover signals as well 

30 as the cross-correlation (CCF) between them. We use the two auto-correlations to deduce the likely shape of the cross- 
correlation, and then fit that 'template' to the cross-correlation to obtain the best estimate of At. 
[0040] The two auto-correlations, ACF CRU and ACF CBU , together with the cross correlation, CCF, may be obtained 
by any of the standard means (see for example in "Random Data: analysis and measurement procedures" by J. S. 
Bendat and A. G. Piersol, Wiley-lnterscience, 1971); we prefer the methods which involve fast Fourier transforms. The 

35 ACF CBU is a good estimate of the 'intrinsic' auto-correlation of the transmitted signals (negligible multipath effect) whilst 
the ACF CRU is likely to be corrupted by multipath signals, i.e. signals delayed by more than about 30 m equivalent 
path. This is shown for a particular example in Figure 6, where trace A is ACF CBU and B is ACF CRU . We can see that 
multipath has given rise to secondary peaks in trace B which appear on both sides of the main peak since the auto- 
correlation is always symmetrical. The 'template' is constructed by taking the left-hand portion of ACF CBU (i.e. the trace 

40 to the left of the central peak in A corresponding to the negative time axis) and joining it to the right-hand portion of 
ACF CRU (i.e. the trace to.the right of the central peak in B corresponding to the positive time axis). This is shown in 
Figure 6 trace C. The measured CCF is shown by trace D, and we can see that the template C reproduces the major 
features. The estimate of the time offset, At, is then obtained by fitting the template to the measured CCF using any 
standard method, allowing the amplitudes of the two halves of the template to be free parameters as well as At. The 

45 fitting procedure also produces a 'goodness of fit' parameter which may be used as a basis for rejection or inclusion 
of the data. 

[0041] In cases where more than three GSM transmitters can be received at both base and rover, it may be an 
advantage to use them all in the position determining process, especially if one or more of the paths is substantially 
indirect. In such a case, the goodness of fit parameter can provide a weighting factor for each Afto be used in evaluating 
so the position. 

[0042] GSM handsets incorporate multi-level Analogue-to-Digital Converters (ADCs) and use adaptive processing 
techniques in order to reduce the errors resulting from multipath propagation. It may therefore be possible to use the 
demodulated signals themselves for the cross correlation function rather than the raw / and Q data streams from the 
receivers, depending on the design of the handset. Such extra delay put in by the adaptive processing must be removed 
55 of course, and no amount of processing can allow for the unknown extra path length of a wholly-indirect signal, but the 
demodulated output may nevertheless be better for estimating time-delays than the pre-processed / and Q raw data 
streams. Shown in dotted outline in Figure 5 is an FM discriminator 19 (FMD). FM discrimination translates a frequency 
offset into a steady (DC) shift in the output, so quadrature signals (I and Q) are not needed, but are replaced by a 



9 



EP 1 118 871 A2 



single signal (J in Figure 5). This offers two further advantages in that (a) half as much data need be transferred over 
the links, and (b) it is unnecessary to adjust for the frequency offset at the position determining processor thus reducing 
the processor overhead. 

[0043] As noted above, several different configurations of a position determining system according to the present 
5 invention are possible, depending on the relative locations of the GSM transmitters (often called Base Transceiver 
Stations, BTSs), the roving receivers (CRUs), the base receivers (CBUs), the position determining processor (CPP), 
and a service provider database (SPD). Three configurations are illustrated in Figures 7a, 7b, and 7c, denoted A, B, 
and C. Each of these is further subdivided into two. Figure 7a shows configuration A in which the CPP and SPD are 
co-located at a position remote from both the CBU and CRU. This could be the case, for example, when one CPP/ 
10 SPD combination serves an entire region. In A1 each BTS has a CBU co-located with it. In A2 the CBU serving the 
local cell embraced by the three BTS is shown somewhere within the cell. Configurations B1 and B2 (Figure 7b) are 
similar to A1 and A2 respectively except that there are now many CPPs. In B1 each BTS has a CBU attached to it; in 
B2 each triangle of three BTSs has a CPP within it of which just one is shown in the Figure. The regional SPD remains 
at a remote location. Configurations C1 and C2 (Figure 7c) are similar to A1 and A2 but with the position processing 
is carried out within each handset. Other configurations are also possible, each having relative advantages making them 
suited to particular commercial applications. 

[0044] In some commercial applications, it may be an advantage to display map information in graphical form as 
part of the response of the service provider database to a request for a position- related service. An example of this is 
given in Figure 8 where it is shown how a map of the local area might be displayed on a screen attached to the position 

20 determining processor or at a monitoring station. The user has requested directions to a location marked by a filled 
circle labelled Z, and the position determining processor has computed the user's position to be somewhere within the 
filled circle labelled W. The user's handset includes a liquid crystal dot-matrix display which can display small portions 
of the larger map as well as character-based information. The first screen displayed on the handset might be as shown 
in the panel marked (a) in Figure 8. This is a large scale map of the local area with the circle, radius 100 m, centred 

25 on the calculated position. The roads are marked 1 , 2, 3 etc., corresponding to the names indicated on the second 
screen, panel (b). The user can scroll back and forth between screens at will. Panel (c) shows a larger scale map of 
the local area, together with a pointer, K, which indicates the direction of the target destination, Z. Clearly, many other 
screens of information are possible, including one giving the address of the target destination, local traffic conditions, 
approximate distance, etc. In a commercial application, there might also be scope for a limited amount of advertising 

30 on the handset display. 

Tests of prototype system 

[0045] A prototype test apparatus using a pair of standard hand portable units and some additional recording equip- 
35 ment was set up in Cambridge, UK. The base station (CBU) was in a three-storey building near the city centre (see 
Figure 9) with the antenna fixed on the roof. The roving receiver (CRU) was carried in a car to various locations in the 
south of the city as shown by the crosses in Figure 9. Three GSM transmitters (BTSs) were used in the trial, located 
respectively in the city centre, at Great Shelford to the south, and at Fulbourn to the south east. Table 2 gives the mean 
positional errors in both co-ordinates, i.e. the differences between the averages of about 10 independent position 
40 determinations and the true positions as measured on a map. Simple data analysis techniques were used with no 
special discrimination against the effects of multipath propagation. The inherent accuracy of the apparatus, found from 
tests with both base receiver and roving receiver side-by-side on the bench, was less than 20 m rms. Tests numbers 
1 to 6 were performed with the rover stationary; tests 7-11 were made at speeds between 40 and 64 km h _1 . 

45 Table 2: 



Location number 


Mean error East (m) 


Mean error North (m) 


1 


66.5 


-59.5 


2 


-5.5 


80.0 


3 


-68.0 


-15.0 


4 


89.0 


90.0 


5 


-19.0 


41.0 


6 


64.0 


79.0 


7 


-73.0 


6.0 


8 


-62.0 


59.0 


9 


17.0 


128.0 


10 


43.0 


31.0 
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Table 2: (continued) 



Location number 


Mean error East (m) 


Mean error North (m) 


11 


95.0 


105.0 



Claims 

1 . A method of estimating the offset (Ar) in time of receipt of a broadcast signal received at two locations (CBU ,CRU), 
10 where the signal received at the one location may have been subject to corruption by multipath effects, the method 

comprising auto-correlating (ACF CBU ), the signals received at the one location (CBU); auto -correlating (ACF CRU ) 
the signals received at the other location (CRU); cross-correlating (CCF) the signals received at the one and other 
locations; constructing a template (trace C Figure 6) comprising that portion of the auto-correlation of the signals 
received at the other location corresponding to the negative time axis, and that portion of the auto-correlation of 
15 the signals received at the one location corresponding to the positive time axis; and measuring the offset (Af) at 

which the template best fits the measured cross-correlation of the signals received at the one and other locations 
as an estimate of the time offset between the signals received at the two locations. 
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